Dynamics and stability of Bose-Einstein solitons in tilted optical lattices 
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Bloch oscillations of Bose-Einstein condensates realize sensitive matter-wave interferometers. We 
investigate the dynamics and stability of bright-soliton wave packets in one-dimensional tilted optical 
lattices with a modulated mean- field interaction g(t). By means of a time-reversal argument, we 
prove the stability of Bloch oscillations of breathing solitons that would be quasistatically unstable. 
Floquet theory shows that these breathing solitons can be more stable against certain experimental 
perturbations than rigid solitons or even noninteracting wave packets. 



PACS numbers: 03.75. Lm; 52.35.Mw; 37.10.Jk 

Matter-wave interferometers are established as stan- 
dard tools for precision measurements of small forces. 
The tiny wavelength of ultracold atoms is an asset, and 
using Bose-Einstein condensates (BECs) can greatly en- 
hance signal-to-noise ratios [1, 2]. In this context, Bloch 
oscillations (BOs) of wave packets in tilted periodic po- 
tentials have been proven very useful [3-6]. Much after 
their early prediction by Zener [7, 8] , BOs were observed 
with electrons in semiconductor superlattices [9, 10], with 
cold atoms in optical lattices [11, 12], and with photons 
in waveguide arrays [13, 14]. BOs are very sensitive to 
dephasing since they rely on coherent Bragg scattering of 
k- vectors from one boundary of the Brillouin zone to the 
other. The slightest lattice imperfection or interaction 
causes random scattering of different fc-components of a 
wave packet, thus broadening its momentum distribution 
and destroying coherent oscillations in real space. 

Long-living BOs of up to 10 cycles with period Tb [15] 
were achieved in a BEC experiment by tuning the scat- 
tering length between Cs atoms to zero with the help of 
a suitable Feshbach resonance [16, 17]. At finite inter- 
action, it appears recommendable to use stable localized 
wave packets, namely, soliton solutions to the nonlinear 
Schrodinger equation, in order to minimize the detrimen- 
tal effects of interaction. Bright solitons arise from a dis- 
persion that counteracts the effect of the nonlinearity. As 
a rule, a soliton can only be stable if the effective mass 
and interaction parameter have opposite signs. In the 
free-space case of positive mass, bright solitons are re- 
alized with atoms that attract each other [18-20]. In a 
lattice, the mass becomes negative close to the edge of 
the Brillouin zone, and there solitons can be prepared 
with repulsive interaction [21]. During BO cycles, the 
mass m(t) changes its sign periodically. Thus, the inter- 
action parameter g(t) has to change accordingly in order 
to respect the criterion of opposite signs; both spatial 
and temporal nonlinearity management schemes to this 
effect have been proposed within the framework of the 
Gross-Pitaevskii equation [22, 23]. 

Independently, we have investigated the stability of 
BOs with an interaction parameter g(t) that is modu- 
lated harmonically in time, and we have found a whole 



family of cases that yield stable BOs of breathing wave 
packets [24]. Conversely, many other instances of g(t) 
result in rapid destruction of BOs. Remarkably, the 
stability criterion developed in Ref. [24] is obviously in- 
compatible with the simple sign rule mentioned above. 
For instance, the Bloch-periodic modulations g(t) = 
±g cos(2-7rf/TB) with opposite signs lead both to stable 
BOs, and this for the same mass m(t). This prompts a 
question that is as fundamentally interesting as it is im- 
portant for practical applications: Is solitonic stability 
helpful to sustain BOs in general? 

In this work we show that stability of BOs is not con- 
ditioned on wave-packet rigidity. We first introduce a 
bounded time in the equation of motion that predicts per- 
fectly stable breathing solutions, in full agreement with 
previous results [ ]. But even under these premises, 
rigid solitons could be expected to be more robust against 
experimental imperfections. We study in detail the rele- 
vant case of a BO decay caused by magnetic fields that 
oscillate off phase. Contrary to expectations, we find 
that the breathing wave packet is more stable than the 
rigid soliton. For the experimentally realistic parameters 
chosen, the rigid soliton is very close to noninteracting. 
Thus, a harmonic modulation of finite interaction turns 
out to effectively stabilize BOs. 

In the mean-field regime, the BEC is described by the 
complex order parameter ^(x,t) in a ID optical lattice 
potential V{x) — V cos 2 (nx/b) with spacing b. If the 
lattice is sufficiently deep, one may use a tight-binding 
approximation, where the condensate is represented by a 
single complex number ^ n (t) at each lattice site [25]. Our 
starting point is thus the nonlinear equation of motion 

iM n = -J(^ n+1 +'$ n - 1 ) + Fbnt> n + g(t)\V n \ 2 t> n (1) 

for the order parameter with normalization = 1- 

Nearest-neighbor sites are coupled by the tunneling el- 
ement J/E w A{V /Ef' 4 exp(-2^V o /.E)/0r", where 
E = h 2 n 2 /2mb 2 is the lattice recoil energy [26]. The 
interaction parameter g = U 2uj x u)yLo z hm /nNa/J is de- 
rived from the scattering length a, the atom number N, 
and the local ground state determined by the lattice fre- 
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quencies u>i [25]. The dispersion relation of this single- 
band model reads e(p) = —2Jcos(pb/h). Its curvature 
or inverse mass mT 1 = 2Jb 2 cos(pb / h) / ti 2 determines the 
wave-packet dynamics under the influence of a constant 
force F that stems from a uniform acceleration of the 
BEC by, e.g., gravity. At fixed lattice geometry, the in- 
teraction parameter g(t) can be controlled by external 
magnetic fields using suitable Feshbach resonances. 

In the following, we use J and b as units of energy and 
length, respectively, and set H = 1. We tackle Eq. (1) by 
separating the rapidly varying Bloch phase p(t)n from a 
smooth envelope A(z, t) comoving with the center of mass 
x{t): # n (i) = e ip ^ n A(n - x(t),t)e^*>. With p(t) = 
-Ft, x(t) = x + 2cos(Ft)/F, and 4> = c/) Q + 2sm(Ft)/F, 
the envelope is found to obey the equation 
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neglecting higher spatial derivatives of A. Note that 
choosing an immobile wave packet with p(0) = as ini- 
tial condition, fixes the phase for the subsequent BOs. 
The inverse mass to(£) -1 = 2cos(Fi) oscillates rapidly, 
but if the interaction is tuned such that m(t)g(t) < 
varies slowly enough, Eq. (2) admits a soliton solution 
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whose quasistatic width is £ = — 2/[g(t)m(t)] > 0. To be 
stable, the soliton must be able to follow this width adi- 
abatically. Otherwise its breathing mode will be driven, 
and other excitations may be created. Therefore, the 
least disruptive way of accelerating a soliton of width 
£o is to impose a perfectly rigid envelope by choosing 
g T (t) = -2m(i)-7£ = g r cos(Ft) with g T = -4/& < 0. 
More extensive studies based on this idea have been put 
forward in Refs. [22, 23]. 

Let us contrast this reasoning with an analytical time- 
reversal argument that we developed after studying the 
stability of BOs under harmonic modulations of g [24]. 
Quite generally, a rigid wave packet is by no means nec- 
essary for persistent BOs. Already for the case g(t) = 0, 
breathing is the rule: In the first quarter of the Bloch 
cycle, the mass is positive and the wave packet spreads. 
When the mass changes sign, the time evolution is re- 
versed and the wave packet recovers its original shape 
at the edge of the Brillouin zone. Thus, the wave packet 
shows periodic breathing on top of the BO, independently 
of its initial shape. Also in the interacting case, one can 
find nontrivial functions g(t) in Eq. (2) compatible with 
this time reversal. Consider the class of periodic func- 
tions 



g(t) = cos(Ft)P(sm{Ft)/F) 



(4) 



where a factor cos(Ft) can be separated from a polyno- 
mial P(r)) in the bounded time variable 
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Time evolution scheme of stable BOs. 
, the interaction parameter of the rigid 



FIG. 1: Left panel: 
The inverse mass m 
soliton <? r as well as that of a breathing soliton gh axe shown 
together with the bounded time g — sin(Ft)/F as function 
of time or momentum — p — Ft. Right panel: The fc-space 
density [obtained by numerical integration of Eq. (1)] of a 
breathing wave packet is a function of g and thus strictly 
periodic in t: the points in time a and b as well as c and d 
show the same distribution, respectively. 



Because dtrj(t) = [2m(i)] _1 , the explicit time dependence 
of the mass factorizes from all terms in the equation of 
motion (2) for A(z,t) = A(z,r](t)): 

id v A(z, rj) = - 3 2 z A(z, V ) + P{g)\A\ 2 A{z, r,). (6) 

The ensuing dynamics for A{z, rf) as function of 77 may 
be quite complicated. However, as rj(t) itself is a peri- 
odic function of time, the solution A(z, t) must also be 
periodic: Any dynamics taking place in the first quar- 
ter of the Bloch period, while 77 runs from to 1/F, is 
exactly reversed in the next quarter, when r\ runs back. 
Figure 1 illustrates this argument by showing the time 
dependence of several key quantities over one Bloch cy- 
cle, as well as a fc-space density plot with clearly visible 
breathing dynamics [30]. 

We stress that Eq. (4) includes both cases g(t) = 
±goCOs(Ft). Although the +cos case does not fulfill 
the soliton stability criterion m(t)g(t) < 0, the preced- 
ing time-reversal argument ensures that both cases lead 
to undamped Bloch oscillations — at least within the ap- 
proximations underlying the equation of motion (2). As 
shown in the right-hand panel of Fig. 1, this prediction is 
confirmed by numerical integration of the tight-binding 
model (1) with a standard fourth-order Runge-Kutta 
method. 

Which of these stable solutions are the most robust 
under variations of experimental control parameters? In- 
deed, even cold-atom experiments suffer from slight im- 
perfections, such as residual uncertainties in the magnetic 
field controlling the interaction term g{t). For instance, 
in the Innsbruck experiment [15] the magnetic field is 
controlled up to lmG. The slope of 61a /G at the zero 
of the Feshbach resonance turns this into an uncertainty 
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FIG. 2: Momentum width (Afc) 2 of a wave packet, Eq. (3), 
with initial spatial width £ = 66.16 (plus a small seed 
noise of 1CP 3 mimicking experimental inhomogeneities), 
Bloch-oscillating in a tilted lattice with F — 2-k/Tb = 
0.15. The interaction parameter is modulated as g(t) = 
goCOs(Ft) + u(t). (a) Random perturbation with frequen- 
cies fi n = nF/N below and above the Bloch frequency 

v(t) = 5^^ =1 Re[a n exp(if2„t)], N = 5. The a n are complex 
random numbers with — 0, |a ?l | 2 = 0.5. (b) Sine pertur- 
bation u(t) = 0.5sin(J ? t). In all cases the lifetime increases 
with the perturbation amplitude go, from the antibreathing 
wave packet go = —10, over the rigid soliton go — g T — —0.06 
and the linear wave packet go — to the breathing soliton 
go = +10. 



Aa = 0.06ao in the scattering length, which is converted 
to the uncertainty of the dimensionless tight-binding in- 
teraction Ag R3 0.4. Note that this uncertainty Ag is 
larger than the interaction g r — — 4/£o needed to create 
a rigid soliton of only moderate width £o > 10. From this 
point of view, realizing a wide rigid soliton is practically 
equivalent to switching the interaction off altogether. 

We study numerically the effect of perturbations of 
g(t) by numerical integration of Eq. (1). Rather than 
the strong force F w 34 in a vertical lattice [15], we 
choose a smaller force F = 0.15, corresponding to a 
slighter tilt. This results in a longer Bloch period and 
a higher sensitivity to dephasing. The broadening of the 
momentum distribution is directly measurable from the 
experimental time-of-flight images and signals the de- 
cay of the wave packet and destruction of BOs in real 
space. In Fig. 2 the fc-space variance (Afc) 2 is shown for 
g(t) = go cos (Ft) + v(t) with two different types of per- 
turbations v(t) and different modulation strengths go- In 
all cases, the momentum distribution starts to broaden at 
a certain time. At a given perturbation, the rigid soliton 
and the linear wave packet, show greater resilience than 
the strongly antibreathing wave packet, but surprisingly 
the breathing wave packet survives even longer. Obvi- 
ously, the +cos modulation stabilizes the Bloch oscilla- 
tions against uncontrolled variations of the interaction 
parameter g. 

Comparing a random superposition of different fre- 



quencies [Fig. 2(a)] to a Bloch-periodic off-phase per- 
turbation proportional to sm(Ft) [Fig. 2(b)], we trace 
back the differences in lifetime to different sensitivities to 
the sine perturbation. In the following, we thus consider 
g(t) = go cos(Ft) + g\ sin(Fi) with an off-phase pertur- 
bation with amplitude g\ of order Ag. In order to un- 
derstand quantitatively why the +cos-modulated wave 
packet can be more robust than the rigid soliton, we 
perform a homogeneous stability analysis. The sudden 
growth of the momentum variance suggests an instabil- 
ity due to the creation of small fluctuations. If these 
perturbations occur on a length scale much shorter than 
the size of the wave packet, the wave packet can be taken 
to be locally homogeneous, |A(z,i)| 2 s» n = l/2£. The 
real and imaginary parts of the small fluctuations to lin- 
ear order, s and d, respectively, obey equations of motion 
that decouple in Fourier modes [24]: 



d\ = -[k 2 cos(Ft) + 2n g(t)]s k , 
Sfe = k 2 cos(Ft)dk- 



(7) 



Thanks to their linearity and time periodicity, these 
equations allow us to apply Floquet theory [28]. In- 
tegrating Eq. (7) over a single period determines the 
Lyapunov exponent Afc that characterizes the exponen- 
tial growth of mode k. For the present case g(t) — 
g cos(Ft) + gisin(Ft), Eq. (7) can be solved analyti- 
cally for gi = 0, using a Bogoliubov transformation 7^ = 
^u>k/k 2 Sk + iy/k 2 /u} k d k , with uj k = yjk 2 (k 2 + 2n g ). 
Perturbation theory to first order in 171 then predicts the 
growth of mode 7^ with 



Jx(2oj k /F) 



(8) 



where J\ is the Bessel function of the first kind. For the 
two cases of primary interest — the rigid and the breath- 
ing soliton — the upper panel of Fig. 3 shows this Lya- 
punov exponent, indistinguishable from the value ob- 
tained by numerical solution of Eq. (7), as function of 
k. Mainly the prefactor k 2 /\u k \ = |1 + 2g n /k 2 \~ 1 
makes the Lyapunov exponents of the breathing wave 
packet (go — 10) smaller than those of the rigid soliton 
(<7o = <7r = -0.06). [31] 

The Lyapunov exponent (8) provides a rather faith- 
ful portrait of the fc-space evolution as obtained by the 
numerics, plotted in the lower panels of Fig. 3. No- 
tably, excitations grow exclusively in the intervals with 
the largest Lyapunov exponents, found at fc-values such 
that k 2 3> n-ogo- Thus, the most unstable mode (indi- 
cated by the vertical line) is found close to fc* defined by 
J[{2k 2 jF) = 0, which gives fc* w 0.96VF « 0.37. The 
predicted growth rate agrees very well with the numerical 
data (inset in the upper panel). These predictions remain 
valid also for smaller go- In the limit go — > 0, the differ- 
ences between different wave packets disappear, and we 
recover the case of a pure sin(i ? t)-modulation that was 
analyzed in [24], and for which Eq. (8) now provides an 
analytical expression. 
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FIG. 3: Decay of a rigid (left) and a breathing (right) soliton 
under the harmonic perturbation of Fig. 2 (b). Upper panel: 
linear stability growth rate, Eq. (8). Lower panel: Strobo- 
scope plot of the fe-space density on a logarithmic color scale. 
Inset: Location and growth of the most unstable mode agree 
with the analytical prediction. 

In essence, this linear-stability analysis applies when- 
ever the excitations are well decoupled in /c-space from 



the original wave packet (the central peak around k = 
visible in Fig. 3). Note that other types of perturbations, 
e.g., a constant offset g 2 as in g(t) = g cos(Ft) + g 2 [29], 
can cause a more homogeneous broadening of the /c-space 
distribution. Its effect is better captured by an ansatz in 
terms of collective variables (cf. [24]). A detailed analy- 
sis is beyond the scope of the present work and will be 
addressed in a forthcoming publication. 

In conclusion, we have connected the physics of soli- 
tons in lattices with a stability analysis of BOs under har- 
monic variations of the interaction. Stability of BOs does 
not rely on soliton stability. Instead, the wave packet 
may start to fall apart, but comes back by virtue of peri- 
odic time reversal. In the presence of instability-inducing 
perturbations, a modulation of the interaction can make 
the wave packet more robust. We explain this behav- 
ior quantitatively via linear stability analysis within Flo- 
quet theory. Finally, let us stress that these results play 
an important role in the design of accurate and reliable 
matter-wave interferometers based on BOs. 
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Eq. (4) covers all stable cases that are Bloch periodic, 

1. e. V2 = 1 in Eq. (8) of [21]. The cases 1/2 / 1 are also 
covered by generalizing the bounded-time argument [27]. 
For go < 0, very small k- values such that k 2 < 2no|(?o| 
have imaginary Bogoliubov frequencies uik, but Eq. (8) 
allows for a smooth analytic continuation, cf. Fig. 3. 



